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Abstract 

Scalar Tensor Vector Gravity(STVG) is one of modified gravity theories developed 
by John Moffat(2005). MOG is abbreviated name for this theory.It can explain a 
galactic rotation curve and the structure formation without dark matter. It can also 
explain acceleration universe without dark energy. But, they obtaion only a spherically 
symmetric, static vacuum solution about MOG. On this theory,the gravitational field 
produced by two point sources is not simply the sum of their respective spherically 
symmetric static vacuum solutions. However,in arXiv:0805.4774, the method to 
adapt MOG to extended distribution of matter is described by phenomenalism. A 
new parameter ^ is introduced in this phenomenalical description. This paper shows 
estimation of MOG's ^ parameter. In coclusion,^ should be less than 0(10^ ) to 
reproduce "flat" rotation curves observed. 



1 Alternative to "Dark side of universe" -Introduction 

The staudard model of cosmology today, the ACDM model, provides an excellent fit to cosmological 
observations. But, most of the composition of the universe is invisible and undetectable which is named 
dark energy or dark matter. This fact provides a strong incentive to seek alternative explanations that 
can explain cosmological observations without dark matter or dark energy. Modified Gravity (MOG) 
[1] has been used successfully to account for galaxy cluster masses [8], the rotation curves of galaxies 
[5][G], velocity dispersions of satellite galaxies, and globular clusters [11]. It can explain the observation 
of the Bullet Cluster [6] without cold dark matter. Besides, MOG also meets the challenge posed by 
cosmological observations. In the paper(arXiv:0710.0364)[4],it is demonstrated that MOG produces an 
acoustic power spectrum, a galaxy matter power spectrum, and a luminosity-distance relationship that 
are in good agreement with observations. 

2 Theory of MOG - Precedent study 
2.1 Action 

The action of Moffat theory is constructed as follows [1]. STVG is formulated using the action principle. 
In the following discussion, a metric signature of [+, — , — , — ] will be used; the speed of light is set to c = 1, 
and we are using the following definition for the Ricci tensor: i?^^ = ^^F^^ - 9^r^„ + T^^T^^^ - ^""^p^lu- 
We denotes the Einstein-Hilbert Lagrangiani/^G = "j^fa + 2A) y/—g, where R is the trace of the 



Ricci tensor, G is the gravitational constant, g is the determinant of the metric tensor g^^^, while A is 
the cosmological constant. Introducing the Proca action Maxwell-Proca Lagrangian for the STVG vector 
field 6,,:. 



La, = --—i^ 



(1) 



where B,, 



da 



dy(j)^, fi is the mass of the vector field, ui determines the strength of the coupling 



between the fifth force and matter, and is a self-interaction potential. The three constants of the 
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theory, G, /.i and w, are promoted to scalar fields by introducing associated kinetic and potential terms 
in the Lagrangian density: 
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where denotes covariant differentiation with respect to the metric 5^^, while Vg, V^, and T4; are the 
self-interaction potentials associated with the scalar fields. The STVG action integral takes the form 
S = J {Cg + £0 + /3s + 'Ca/) d'^x, where Cm is the ordinary matter Lagrangian density. 



2.2 Weak field approximation of MOG 

The field equations of STVG can be developed from the action integral using the variational principle. 
First a test particle Lagrangian is postulated in the form/Ixp = ^™ + c«uj(j5^^u^, where m is the test 
particle mass, a is a factor representing the nonlinearity of the theory, is the test particle's fifth-force 
charge, and = dx^ /ds is its four-velocity. Assuming that the fifth-force charge is proportional to 
mass,q5 = Km, the value of k = -^/Gtv/o; is determined and the following equation of motion is obtained 
in the spherically symmetric, static and weak gravitational field of a point mass of mass M : 



GnM r 



l + a-a(l + r/A)e-''/^ 



(3) 



where Gn is Newton's constant of gravitation. Further study of the field equations allows a determination 
of a and A for a point gravitational source of mass M in the form: 
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^ fg^-lV (5) 
{^fM + Ef \Gn J 

1/2 — 

The constants D and E various astronomical observation yield the following values: D ~ 6250Mq' kpc 

1/2 

E ~ 25000AfQ , In the weak-field approximation, STVG produces a Yukawa-like modification of the grav- 
itational force due to a point source. Nearby a source gravity,the repulsive force which comes from vector 
field counteracts attraction. But far from a source gravity, Yukawa-like repulsive force cannot reach be- 
cause vector field is massive. Intuitively, this result can be described as follows: far from a source gravity 
is stronger than the Newtonian prediction, but at shorter distances, gravity is comparable to Newtonian. 



2.3 Adaption MOG to extended distribution of matter 

The gravity acceleration equation of the STVG in a weak gravitational field is shown in the Eq.3. But, This 
equtaion is adaped only gravitational field of "a point mass" . It is not applied simple superposition 
principle of forces without self-contradiction. In a theory that offers linear behavior in the weak field 
limit, for example general relativity ,an extended distribution of matter can be well approximated by a 
large number of point masses. However MOG is not so: the gravitational field produced by two point 
sources is not simply the sum of their respective spherically symmetric static vacuum solutions. Presently, 
wc do not have solutions to the MOG field equations for extended matter distributions. But, they slove 
this problem phenomenologically in in Moffat ct al.(2009)[ I •'!]. They seek an effective mass function 
Mef f{x,r), to be used in place of M in (5) and (4), that determines an "effective mass", to be used in 
the formulae for a and A. The function should yield the mass of the source in the case of a point source, 
and a mass proportional to volume in the case of a constant distribution. One function that satisfies 
these criteria is in the form, 

M,ff{x,r) = I pCx)cxp d% (6) 

where ^ is a parameter which should be decided to satisfy a astronomical phenomenon observed. 
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3 Estimation of ^ parameter by N-body calculation - Method 

We let the Eq.6 become disintegration for N-body simulation. And the following expressions are derived by 
adapting itself to many body system. In the following expression, "i-particle" is gravity source particle, "j- 
particle" is received particle. 



GeffiiJ) 
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(8) 



This calculation needs the information of all the particles whenever we calculate force between two 
particles. That is why the calculation number of times becomes 0{N^). By the normal method, we 
cannot carry out such an enormous calculations. But, We use "mass shell and G shell aproximation" in 
our precedent study [IG]. Then, we can reduce the calculation number of times to 0{N'^). To estimate 
value of ^ which can reproduce observational phenomena ,we calculated the acceleration which each 
particle of the expnential disk. The following is parameter of the disk.Thease are comparable to the 
Milky way. This ^ is the non-dimensional parameter which occured from phenomenological method. 



IR(radius) 25kpc 

IM(mass) 5 x 10^°Mq 

IT(time) 0.26Gry 

IV(velocity) 92km/s 

Table 1: The standardization of the unit 



i?d (scale length of disk) 


2.5kpc 


i?cut (radius of disk) 


25.0kpc 


(scale length of disk height) 


O.Skpc 


Md (total mass of disk) 


5 X 10^°Mq 


eps (softening length) 


0.048kpc 


Nbody(number of particles) 


10000 



Table 2: Parameters of the disk 



Theoretically,it is no problem that either ^ = 10 or ^ = 10^".Then,we can't guess the rough size of ^. 
We perform order-estimation about We examined 7 models:^ = 10^^, 10^^, 10°, 10^, 10^, 10^, 10^. 
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Figure 1: The ratio of central acceleration by MOG and the central acceleration by the Newton gravity. We 
denote this as effective gravitational constant. 
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Figure 2: The prediction of the galaxy rotation curves under dynamical equilibrium. 
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4 ^ is less than 0(10^) - Result and consideration 

TheUjWe calculated the MOG acceleration of each particle in this disk. The result is shown in Fig.l. This 
figure shows the ratio of central acceleration by the modified gravity and the central acceleration by the 
Newton gravity as an effective gravitational constant. ^ in less than O(10^),as much as ^ grows large, the 
effective gravitational constant in galaxy disk scale becomes strong. If ^ is larger than C'(fO'^),as much 
as ^ grows large, the effective gravitational constant in galaxy disk scale becomes week. In the case of 
^ = 10^, the effective gravitational constant accords with Newton gravity. 

Fig. 2 shows the prediction of the galaxy rotation curves under dynamical equilibrium about each 
model. For a comparison, we calculated the prediction of the galaxy rotation curves on Newton gravity 
about the same disk. The result whose ^ is larger than O{l0^) make a little difference than Newton 
prediction. 

Considering the flat rotation curve observed,^ should be less than O(IO^). 

5 Future work 

Now, we are carrying out a large-scale calculation about galactic dynamics and structure formation of 
universe. From a dynamical viewpoint, we may obtain a further limitation of ^. Nurmerical calculation 
are carried out on Cray XT4 at Center for Computational Astrophysics(CfCA) of National Astronomical 
Observatory of Japan and SRI6000 at YITP in Kyoto University. 
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